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Abstract 

Enviroment - caused dissipation disrupts the hamiltonian evolution of all quantum systems not 
fully isolated from any bath. We propose and examine a feedback-control scheme to eliminate 
such dissipation, by tracking the free hamiltonian evolution. We determine a driving-field that 
maximizes the projection of the actual molecular system onto the freely propagated one. The 
evolution of a model two level system in a dephasing bath is followed, and the driving field that 
overcomes the decoherence is calculated. An implementation of the scheme in the laboratory using 
feedback control is suggested. 
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The Schrodinger equation describes time evolution of isolated quantum systems. Am 
real system, unless very carefully isolated, will decohere due to bath interactions^, ODD 

□□□□□□ 



originated from the idea of changing the local field in order to follow a 
fl fi Q Q. TVacMn g can b e considered as local control with a time 



The approach to decoherence correction taken here derives from tracking control theory. 
Tracking contro 
reference state 

dependent target. Since generation of interference pathways is a global process in time, an 
intervention at one instant influences the subsequent evolution and can interfere at later 
times with another intervention. 

The basic idea here is to use a freely propagating molecule as a tracking target for an 
identical molecule in a dissipative environment. The perfect tracking field will then null the 
dissipative forces and protect the molecule from decoherence. The question to be addressed is 
how can an external field, that generates a unitary transformation, protect the system against 
non unitary dissipative forces? A computational model of an on-line, dynamic, decoherence 
control method is explored. The method is established by simultaneously following the 
evolution in time of a system with and without the bath, thus using the overlap between 
the two systems to construct a correcting electric field which is directly applied back on the 
system. The method is explored for the fastest dissipative mechanism, environment-induced 
electronic quenching from an excited state to a lower state (the approximate decay time in 
solid state systems is less than 50 fsec). The coherence time of the system is extended by 
an order of magnitude using the simplest tracking method. 

Measuring the overlap between the controlled and the target system also underlies a 
proposed experimental approach. This observable can then be employed as a basis for a 
learning loop used to optimize a control field. 

The molecular system is described by the density operator p, a function of nuclear and 
electronic degrees of freedom. The reference or target system p tar , evolves freely under its 
Hamiltonian Ho: 

'Ptar = -«[H , Ptar] (!) 

where H = 1 is chosen. The subjected system p s is coupled to a bath and therefore evolves 
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according to the dynamics of an open quantum system |l 

ps = -i[n ,p s ] +£d{Ps) 



(2) 



Here describes the bath-induced decoherence dynamics, conveniently cast into the 
Lindblad semigroup form: 

£d(Ps) = £ (f^sF] - P S }) , (3) 

where {A, B} = AB + BA is the anti commutator and F are operators from the Hilbert 
space of the system. The nature of the bath interaction determines the form of the Lindblad 
operators F. The choice of F determines the dissipative model considered. 

Our purpose is to control the system coupled to the bath and force it to follow as closely 
as possible the dynamics of the freely evolving reference system, p tar . The external control 
field e(t) is coupled to the system through the transition dipole of the molecule. The control 
Hamiltonian becomes, H c = fie(t) leading to the dynamics of the controlled system p c : 

pc = -i[Ho, Pc) ~ «[H C , pc] + C D {pc) ■ ( 4 ) 

The challenge of decoherence control is to find the field, e(t), that induces the dynamics 
of the system p c (t) to be as close as possible to p tar . This is attacked by maximizing the 
overlap functional J(t) between the target and controlled states: 

J(t) = (Pc ■ Ptar) = Tr{p c p tar } = (p tar ) (5) 

J(t) can be interpreted as the expectation of the target density operator in the controlled 
state. To achieve this target the control field should increase the expectation of the target 
operator. The Heisenberg equation of motion for the observable (p tar ) generated by the 
control Hamiltonian becomes: 

J t (Ptar) = -«([H c ,Pter]) = (6) 

-ie(t)Tr{ppt ar pc - PtarAPc} • (7) 

The control field is constructed by requiring maximal J(t), so that 4(Pt ar ) > at any 
instant, leading to: 

e(t) = -iK Tr{pcPPtar - PtarApc}* , (8) 
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(9) 



where K = K(t) is a positive envelope function with dimension • From Eq. (JHJ) it is 

clear that the correcting field e(t) = if pc approaches Ptar- 

We adopt a molecular model system with two electronic states described by the density 
operator: 

P Peg 
Pge P g 

where the indices g and e designate the ground and excited electronic states and the sub 
matrices are functions of the nuclear coordinates. The Hamiltonian of the system has the 
form: 

11,, ( " M • (10) 

with H g / e = T + V g / e . T is the kinetic energy operator, V g and V e are the potential energy 
operators on the ground and excited electronic states and V eg is the nonadiabatic coupling 
potential. The control Hamiltonian is chosen as: 

\-fke(t) J 

where fit, is the coordinate dependent electronic transition dipole element. 

The dissipation operator F q is chosen to model fast electronic quenching as F q = 
y/%( I e) {g I ) . This choice of dissipative term will lead to a process with a characteristic 
quenching rate 7 : 

Cn{p) = l\* e ~ A?eV2 I • (12) 

\ "Peg/2 P e 

The control field form Eq. (jHJ) becomes: 



e(t) = KIm(Tr Q {prfi^ ar } + Tr Q {plfi^ r } 
-Tr Q {pf ar fiLpl} - Tr Q {pt ar ppf} 
-Tr Q {p 9 tar fiip^} + Tr Q {pfiipZ} 
-Tr Q {pflfiipf} + Tr Q {p^fiip e tar } ) 



(13) 



where Trg is a partial trace over the coordinates. 

The computational model is constructed from two electronic states of a diatomic molecule 
represented by two diabatic potential energy surfaces that are approximated by quadratic 
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Nuclear Coordinate 



FIG. 1: The potential energy surface of the molecular model and the ground initial state. The parameters 
describing the model are: uu g = uj e = 0.07 eV, A = .7 eV ,V ge / eg = 0.05 eV, Q 9 = , Q e = —0.1 and 
/} = 1.0. 



functions with linear interstate couplings [16]. Using dimensionless normal coordinates: 

V 9 (Q) = -A + f(Q-Q,) 2 , 

K(Q ) = A + f (Q-Q e ) 2 , 
where u> 9 and c«j e are the vibrational frequencies of the ground and excited electronic states, 
2A is the adiabatic excitation energy. Q e and Q 9 are the equilibrium bond lengths of the 
two electronic states. 

The time evolution was obtained by solving the time dependent Liouville-von Neumann 

equation. A grid was used for the spatial coordinates and time propagation used the Cheby- 

chev scheme Il7|. The initial ground vibronic state of the system was calculated via a 
I I 

relaxation scheme . A pump pulse transfers a significant fraction of the population from 
the ground to the excited electronic state. A Gaussian form is chosen for this pulse: 



€ 



pump 



(t) = ^-C-^JVKe^' , (15) 



The carrier frequency ujl is chosen to match the difference between the ground and excited 
electronic potentials at the minimum of the ground state Q = 0. The width of the pulse 0£ 
was adjusted to a FWHM duration of 12 fsec and the amplitude eo • [i = 0.228eU. 

At this point three parallel propagations of the density operator were performed. The 
first for p tar was carried out only with the Hamiltonian term Eq. (jj). The system density 
operator p s , was propagated with the dissipative Co included, Eq. 0. The controlled 
density operator, p c was propagated with both the dissipative Cd and the correcting field 
applied terms, Eq. (JU). The correction field was calculated according to Eq. (fT3j) at each 
time step and fed back in the next time step. 
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FIG. 2: The population of the excited state as a function of time for the target state p tar (red), the 
controlled density operator p c (black) and the uncontrolled system density operator p s (green) for different 
magnitudes of quenching(columns) and A's(rows). 7 = O.OOS/sec -1 . 



The correction scheme was tested for the the system shown in Fig. ^ The free dynamics 
represents a complex population oscillation between the two electronic states. The excited 
state population is shown in Fig. |2]for the controlled p c , uncontrolled p s and target p tar for 
different quenching parameters 7 and energy gap A. In all cases a fast decay of the excited 
state is seen for the uncontrolled state. The controlled state is able to track the target 
state and maintain the population. When the quenching increases the controlled system 
is still able to follow the overall dynamics of the target but not the finer details. As can 
be seen in Fig. |2l the control ability increases with the energy gap A. The control field 
e is shown in Fig. |3] in time and frequency. In general the field is composed of a central 
frequency corrosponding to the energy gap modulated by the vibrational frequency. The 
actual vibronic lines are missing from the spectrum [19] . Wigner plots ( not shown ) show 
phase locking between the frequency components. 

The changes in purity Tr{p 2 } and in the scalar product with the target state (p • p tar ), 
show a different viewpoint on the decoherence control Fig. HJ The uncontrolled state p s 
undergoes fast exponential decay of both measures in time. The controlled purity Tr{p 2 J } 
maintains a high value as does as the scalar product (p c ■ p tar ) The general trend is a 
linear decay in time. The oscillations around this d ecay and the temporal increase in purity 



suggest that a cooling mechanism is taking place 



19j. To gain insight on this possibility 
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FIG. 3: Typical control field in frequency (top left) and time domain (bottom) for 7 = 0.003(red) and 
7 = 0.006 /sec" 1 (blue). 



the correcting field was stopped for 100 fsec and resumed again, Fig. 03 As expected, once 
the control field is turned off the purity and the scalar product of the controlled state p c 
decease sharply. When the controlling field is resumed an almost linear increase in both 
parameters can be observed. Such an effect can only be a result of a cooling mechanism 
which represents an interplay between the nonunitary dissipation and the unitary control 



field 



It is clear 



20j that simple unitary transformation cannot induce cooling. 



The present results show that the decoherence can be suppressed, extending good over- 
lap with the target by an order of magnitude in time. The model calculation establishes 
the principle that a correcting external field can be found which protects the system from 
decoherence. Here we discuss control to overcome fast electronic dephasing; while the rapid 
modulations in Fig. El are not currently attainable, the essential behaviour is correct. 

The proof of principle then shifts the attention to finding such a field in the laboratory. 
The idea is to use an actual freely evolving hamiltonian subsystem as the target for an 
identical hamiltonian subsystem subject to a dissipative environment such as in solution or 
on a surface. The correcting feedback should be set such that the difference in transient 
absorption of the two molecules is minimized. 

To implement the present scheme in the laboratory the key element is to apply a scalar 
product between the states of the reference and controlled systems: (p c ■ p tar ) (displayed 
in FigEJ. To make this task possible the two isolated states should become part of the 
same quantum system. This can be carried out by placing the controlled system p c and 
the reference system p tar in two branches of an interferometer. A pair of twin ultrashort 
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FIG. 4: Top - The purity Tr{p 2 } (solid) and the scalar product with the target state (p • p tar ) (dashed) for 
the target state p c (red) the controlled density operator p c (black) and the reference uncontrolled density 
operator p s (green). Bottom - The population of the excited state, same color codes. 
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FIG. 5: Turning off the control field. Top -The purity Tr{p 2 } (solid) and the scalar product with the target 
state (p • p tar ) (dashed) for the target state p tar (red) the controlled density operator p c (black) and the 
reference uncontrolled density operator p s (green). The vertical lines indicate the time interval where the 
control field is turned off. Bottom - The population of the excited state as a function of time, same color 
codes. 



photons which can be created by down conversion are split spatially. Then they are directed 
to the two branches of the interferometer interrogating the two systems. The transmitted 
photons are then redirected to interfere with each other. Any difference between the target 
and the controlled system will degrade this interference. The interference signal can then 
be used as a feedback to generate the correcting field 
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Enviroment - caused dissipation disrupts the hamiltonian evolution of all quantum systems not 
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feedback control is suggested. 
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The Schrodinger equation describes time evolution of isolated quantum systems. Any 
real system, unless very carefully isolated, will decohere due to bath interactions[l-10]. 

The approach to decoherence correction taken here derives from tracking control theory. 
Tracking control originated from the idea of changing the local field in order to follow a 
reference state [11-14]. Tracking can be considered as local control with a time dependent 
target. Since generation of interference pathways is a global process in time, an intervention 
at one instant influences the subsequent evolution and can interfere at later times with 
another intervention. 

The basic idea here is to use a freely propagating molecule as a tracking target for an 
identical molecule in a dissipative environment. The perfect tracking field will then null the 
dissipative forces and protect the molecule from decoherence. The question to be addressed is 
how can an external field, that generates a unitary transformation, protect the system against 
non unitary dissipative forces? A computational model of an on-line, dynamic, decoherence 
control method is explored. The method is established by simultaneously following the 
evolution in time of a system with and without the bath, thus using the overlap between 
the two systems to construct a correcting electric field which is directly applied back on the 
system. The method is explored for the fastest dissipative mechanism, environment-induced 
electronic quenching from an excited state to a lower state (the approximate decay time in 
solid state systems is less than 50 fsec). The coherence time of the system is extended by 
an order of magnitude using the simplest tracking method. 

Measuring the overlap between the controlled and the target system also underlies a 
proposed experimental approach. This observable can then be employed as a basis for a 
learning loop used to optimize a control field. 

The molecular system is described by the density operator p, a function of nuclear and 
electronic degrees of freedom. The reference or target system p tar , evolves freely under its 
Hamiltonian H : 



where h = 1 is chosen. The subjected system p s is coupled to a bath and therefore evolves 
according to the dynamics of an open quantum system [15]: 



Ptar 



i [ H 0, Ptar] 



(1) 



Ps 



i[H ,p s ] +C D (p s ) 



(2) 
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Here C D describes the bath-induced decoherence dynamics, conveniently cast into the 
Lindblad semigroup form: 

£d(Ps) = E (*iPs*i " &}) > ( 3 ) 

where {A, B} = AB + BA is the anti commutator and F are operators from the Hilbert 
space of the system. The nature of the bath interaction determines the form of the Lindblad 
operators F. The choice of F determines the dissipative model considered. 

Our purpose is to control the system coupled to the bath and force it to follow as closely 
as possible the dynamics of the freely evolving reference system, p tar . The external control 
field e(t) is coupled to the system through the transition dipole of the molecule. The control 
Hamiltonian becomes, H c = p,e(t) leading to the dynamics of the controlled system p c : 

p c = -i[H 0) p c ] ~ i[H C) p c ] + C D {p c ) ■ ( 4 ) 

The challenge of decoherence control is to find the field, e(t), that induces the dynamics 
of the system p c (t) to be as close as possible to p tar . This is attacked by maximizing the 
overlap functional J(t) between the target and controlled states: 

J{t) = (Pc ■ Ptar) = Tr{p c p tar } = (p tar ) (5) 

J(t) can be interpreted as the expectation of the target density operator in the controlled 
state. To achieve this target the control field should increase the expectation of the target 
operator. The Heisenberg equation of motion for the observable (p tar ) generated by the 
control Hamiltonian becomes: 

J t (Ptar) = - «([H c ,Ptar]) = (6) 

-ie(t)Tr{p,pt ar pc ~ PtarApc} • (7) 

The control field is constructed by requiring maximal J(t), so that ^(p tar ) > at any 
instant, leading to: 

e(t) = -%K TrjpcPPtar ~ PtarPPc}* , (8) 

where K = K(t) is a positive envelope function with dimension ■ From Eq. (8) it is 

clear that the correcting field e(t) = if pc approaches ptar- 
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We adopt a molecular model system with two electronic states described by the density 
operator: 



(9) 



where the indices g and e designate the ground and excited electronic states and the sub 
matrices are functions of the nuclear coordinates. The Hamiltonian of the system has the 
form: 

Ho= ( He Vey \ , (10) 

\y 9 e H g 

with Hg/ e = T + V g / e . T is the kinetic energy operator, V g and V e are the potential energy 
operators on the ground and excited electronic states and Y eg is the nonadiabatic coupling 
potential. The control Hamiltonian is chosen as: 

tu=( -^ {t) ) , (ID 

where p, is the coordinate dependent electronic transition dipole element. 

The dissipation operator F q is chosen to model fast electronic quenching as F q = 
y/%(\e){g\). This choice of dissipative term will lead to a process with a characteristic 
quenching rate 7 : 

Cd{p) = 7 I ^ ~* 9 ° /2 1 • (12) 



-Peg/ 2 Pe 



The control field form Eq. (8) becomes: 



-TrqipTarP-pVs ~ TrqiplarP-PT} 
-TrQ{pt ar P'pl g } + Tr Q {f% p,p t 9 ar } 
-Tr Q {pl%(ipf} + Tr Q {p^p,p e tar } ) 



(13) 



where Ttq is a partial trace over the coordinates. 

The computational model is constructed from two electronic states of a diatomic molecule 
represented by two diabatic potential energy surfaces that are approximated by quadratic 
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Nuclear Coordinate 

FIG. 1: The potential energy surface of the molecular model and the ground initial state. The parameters 
describing the model are: u g = to e = 0.07 eV, A = .7 eV ,V ge /e g = 0-05 eV, Q g = , Q e = —0.1 and 

A = l.o. 

functions with linear interstate couplings [16]. Using dimensionless normal coordinates: 

V S (Q) = -A + |(Q-Q 9 ) 2 , 

V e (Q ) = A + f (Q-Q e )\ 
where ui g and u e are the vibrational frequencies of the ground and excited electronic states, 
2A is the adiabatic excitation energy. Q e and Q g are the equilibrium bond lengths of the 
two electronic states. 

The time evolution was obtained by solving the time dependent Liouville-von Neumann 
equation. A grid was used for the spatial coordinates and time propagation used the Cheby- 
chev scheme [17]. The initial ground vibronic state of the system was calculated via a 
relaxation scheme [18]. A pump pulse transfers a significant fraction of the population from 
the ground to the excited electronic state. A Gaussian form is chosen for this pulse: 

WW = e e-^— ) 2 /M e ^* , (15) 

The carrier frequency lol is chosen to match the difference between the ground and excited 
electronic potentials at the minimum of the ground state Q = 0. The width of the pulse a L 
was adjusted to a FWHM duration of 12 fsec and the amplitude e ■ /j, = 0.2286V. 

At this point three parallel propagations of the density operator were performed. The 
first for p tar was carried out only with the Hamiltonian term Eq. (1). The system density 
operator p s , was propagated with the dissipative Cp included, Eq. (2). The controlled 
density operator, p c was propagated with both the dissipative Cd and the correcting field 
applied terms, Eq. (4). The correction field was calculated according to Eq. (13) at each 
time step and fed back in the next time step. 
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FIG. 2: The population of the excited state as a function of time for the target state p tar (red), the 
controlled density operator p c (black) and the uncontrolled system density operator p s (green) for different 
magnitudes of quenching(columns) and A's(rows). 7 = 0.003/sec -1 . 

The correction scheme was tested for the the system shown in Fig. 1. The free dynamics 
represents a complex population oscillation between the two electronic states. The excited 
state population is shown in Fig. 2 for the controlled p c , uncontrolled p s and target p tar for 
different quenching parameters 7 and energy gap A. In all cases a fast decay of the excited 
state is seen for the uncontrolled state. The controlled state is able to track the target 
state and maintain the population. When the quenching increases the controlled system 
is still able to follow the overall dynamics of the target but not the finer details. As can 
be seen in Fig. 2, the control ability increases with the energy gap A. The control field 
e is shown in Fig. 3 in time and frequency. In general the field is composed of a central 
frequency corresponding to the energy gap modulated by the vibrational frequency. The 
actual vibronic lines are missing from the spectrum [19]. Wigner plots ( not shown ) show 
phase locking between the frequency components. 

The changes in purity Tr{p 2 } and in the scalar product with the target state (p • p tar ), 
show a different viewpoint on the decoherence control Fig. 4. The uncontrolled state p s 
undergoes fast exponential decay of both measures in time. The controlled purity Tr{p 2 c } 
maintains a high value as does as the scalar product (p c ■ p tar ) The general trend is a 
linear decay in time. The oscillations around this decay and the temporal increase in purity 
suggest that a cooling mechanism is taking place [19]. To gain insight on this possibility 
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